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Degeneracy  in  the  Ensemble  Monte  Carlo  Method  for 
High-Field  Transport  in  Semiconductors 

P.  LUGLI  AND  D.  K.  FERRY,  senior  member,  ieee 


Abslmcl—An  algorithm  to  include  the  Pauli  exclusion  principle  in  the 
Ensemble  Monte  Carlo  method  is  presented.  The  results  indicate  that 
significant  changes  in  the  transport  properties  of  GaAs  have  to  be  ex¬ 
pected  when  degenerate  conditions  are  reached.  Important  repercus- 
siims  should  be  found  in  the  modeling  of  microwave  devices,  where  one 
often  deals  with  highly  doped  regions. 


I.  Introduction 

The  ensemble  monte  carlo  (EMC)  method 

has  been  widely  used,  in  recent  years,  in  connection 
with  the  study  of  nonlinear  transport  in  semiconductors 
[1],  (2).  It  has  also  been  successfully  applied  to  the  sim¬ 
ulation  of  semiconductor  devices  [3],  14].  The  EMC 
method  is  a  semiclassical  technique,  in  that  it  simulates 
electrons  as  classical  panicles  undergoing  scattering  events 
that  are  calculated  according  to  quantum  mechanical  tran¬ 
sition  probabilities.  Classical  statistics  are  used,  and  the 
Pauli  exclusion  principle  is  neglected.  On  the  other  hand, 
many  of  the  situations  of  interest  (for  example,  in  micro- 
wave  devices)  involve  high  concentrations  of  electrons  and 
impurities,  where  degeneracy  is  expected  to  play  an  im¬ 
portant  role.  In  the  following,  we  propose  an  algorithm  to 
include  the  Pauli  exclusion  principle  in  the  EMC  tech¬ 
nique.  Its  effect  on  the  transport  characteristics  of  GaAs 
will  be  discussed.  While  attempts  have  previously  been 
made  to  include  degeneracy  in  a  single-particle  Monte 
Carlo,  these  have  had  only  limited  application  [5].  Here, 
we  show  that  the  unique  properties  of  the  EMC  allow  a 
•Straightforward  inclusion  of  degeneracy. 

The  influence  of  the  exclusion  principle  on  the  average 
transport  properties  of  semiconductors  can  be  of  great  im¬ 
portance.  For  example,  when  modeling  of  microwave  de¬ 
vices  is  treated,  one  often  deals  with  highly  doped  regions 
that  are  expected  to  be  degenerate.  No  Monte  Carlo  sim¬ 
ulation  of  actual  semiconductor  devices  performed  so  far 
has  taken  the  Pauli  principle  into  account  even  for  con¬ 
centrations  well  above  lO''  cm'^  [6].  Furthermore,  when 
balance  equations  are  used  [7],  it  is  always  assumed  that 
electrons  in  the  low  field  regions  (typically  characterized 
by  high  doping  levels)  have  the  same  average  kinetic  en¬ 
ergy  as  the  lattice,  in  terms  of  their  characteristic  tem- 
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perature.  The  results  of  the  degenerate  EMC  calculations 
indicate  that  a  proper  account  of  the  fermion  statistics 
causes  significant  changes,  especially  to  the  electron 
“temperature”  (the  average  energy  of  the  carriers  is  re¬ 
lated  now  to  Ef,  and  not  to  ksT). 

Although  the  EMC  is  a  semiclassical  technique  that 
simulates  electrons  moving  as  classical  particles,  in  reality 
the  electrons  behave  like  fermions.  They  must  still  obey 
the  Pauli  exclusion  principle:  each  quantum  nuipber  is 
available  for  at  most  two  electrons,  which  differ  in  their 
spin  quantum  number.  Following  the  simple  theory  of  the 
Fermi  gas  [8].  the  number  of  available  it -space  states  for 
a  system  enclosed  in  a  volume  V  is  given  by  K/(2t)^.  At 
zero  temperature,  the  electrons  All  up  a  region  of  it -space 
called  the  Fermi  sphere.  The  radius  of  this  sphere,  the 
Fermi  wave- vector  kf,  is  related  to  the  electron  cencen- 
tration  n  (at  absolute  zero)  by 

lit,l  =  (1) 

where  n  is  the  net  electron  concentration.  At  nonzero  tem¬ 
peratures,  electrons  spill  out  of  the  Fermi  sphere,  moving 
into  states  previously  vacant  and  thereby  allowing  for  cur¬ 
rent  flow  to  occur.  In  metals,  only  the  electrons  in  a  region 
of  the  order  of  kg  Tg  (where  Tg  is  the  ambient  temperature) 
from  the  Fermi  level  can  take  part  in  collision  and  con¬ 
duction  processes.  Similar  effects  happen  in  semiconduc¬ 
tors  under  degenerate  conditions,  e.g.,  when  the  Fermi 
level  lies  in  the  conduction  band.  In  GaAs  at  300  K,  this 
occurs  lorn  >  4.6  X  lO'^  cm~'.  The  Pauli  exclusion  prin¬ 
ciple  can  be  thought  of  as  a  many-body  effect  that  influ¬ 
ences  the  transport  properties  of  degenerate  semiconduc¬ 
tors  by  limiting  the  phase  space  available  for  electronic 
transitions. 


II.  Approach 

The  probability  of  an  electronic  transition  from  a  state 
k  into  a  state  k'  is  in  general  proportional  to  the  proba¬ 
bility  f(k)  that^the  initial  state  is  occupied  and  the  prob¬ 
ability  1 1  -/(?')]  that  the  final  s^te  is  unoccupied.  J'hus, 
the^ SMUerijig  rate  from  k  lo  k'  is  th^n  P(k,  k')  = 
S{k,  k')f{k)  [1  -/(?')!,  where  S(k,  k')  is  the  transi¬ 
tion  probability,  usually  calculated  from  the  Fermi  Golden 
Rule.  The  standard  Monte  Carlo  procedure  works  within 
the  approximation f(k')  =  0.  That  is,  it  considers  all  final 
states  as  being  available.  The  latter  method  is  only  appli- 
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cable  to  nondegenerate  conditions,  where  the  effect  of  the 
Pauli  principle  is  negligible.  The  problem  of  working  with 
P(k,  k')  in  the  Monte  Carlo  algorithm  lies  in  the  fact  that 
f(k')  is  not  known  a  priori.  However,  with  the  EMC  tech¬ 
nique,  the  distribution  function  evolves  with  the  ensemble 
and  is  known  at  every  instant  of  time.  One  possibility  that 
allows  us  to  overcome  the  problem  posed  by  the  fector 
|1  -/(?')]  is  to  use  the  differential  scattering  rates  rather 
than  the  integrated  rates  in  the  simulation  (IJ.  P(k.  k') 
can  then  be  calculated  for  discrete  values  of  the  initial  and 
final  momentum  and  updated  at  fixed  time  intervals  dur¬ 
ing  the  simulation.  Storage  requirements,  however,  will 
cause  a  severe  limitation  of  the  efficacy  of  this  approach, 
and  this  suggests  a  different  approach. 

An  alternative  method  makes  use  of  the  procedure  sug¬ 
gested  for  a  one-particle  Monte  Carlo  (S).  To  the  best  of 
our  knowledge,  this  early  work  has  been  the  only  attempt 
to  include  the  Pauli  exclusion  principle  in  Monte  Carlo 
calculation,  and  has  not  become  used  since  the  early  re¬ 
port.  This  method  is  extended  here  to  the  EMC  technique, 
where  it  reaches  real  applicability.  The  self-consistent  al¬ 
gorithm  proposed  by  Bosi  and  Jacoboni  resorts  to  the  re¬ 
jection  technique  |9|  to  account  for  the  occupation  prob¬ 
ability  of  the  final  state  at  each  scattering  event.  At  the 
time  of  the  selection  of  the  particular  scattering  process, 
the  final  state  of  the  transition  is  still  not  known.  The  ex¬ 
clusion  principle  can  therefore  be  neglected  for  the  deter¬ 
mination  of  the  length  of  the  free  flight,  and  for  the  choice 
of  the  scattering  mechanism  and  the  appropriate  final 
state.  However,  once  the  final  state  is  selected, /(A')  is 
known  from  the  EMC  even  during  the  transient  phase  and 
a  random  number  between  0  and  1  can  be  u.sed  to  accept 
or  reject  the  transition.  As  the  procedure  is  iterated,  a 
steady  state  will  eventually  be  reached,  with  the  proba¬ 
bility  of  transitions  into  the  state  k '  proportional  to  the 
occupancy  of  that  state.  In  principle,  the  distribution  func¬ 
tion  is  known  already  with  the  EMC  and  could  be  incor¬ 
porated  into  the  scattering  rates,  as  suggested  above. 
However,  the  continuing  need  to  recalculate  the  scattering 
rates  with  this  latter  approach  makes  the  rejection  tech¬ 
nique  a  much  more  efficient  numerical  technique. 

The  most  delicate  point  of  this  method  involves  the  nor¬ 
malization  of  the  distribution  function/f  k).  Therefore,  the 
extension  of  the  method  to  the  EMC  algorithm  will  be  de¬ 
scribed  in  some  detail.  We  take  N  to  be  the  number  of 
electrons  simulated  in  the  ensemble  and  n  the  real  electron 
concentration.  Then  the  effective  volume  V  of  “real 
space"  that  is  effectively  occupied  by  this  A-particle  sys¬ 
tem  is  V  =  NIn.  The  density  of  allowed  wave-vectors  of 
one  spin  in  A -space  is  given  by  K/(27r)'.  A  grid  is  set  up 
in  the  three-dimensional  wave-vector  space,  with  an  ele¬ 
mentary  cell  whose  volume  is  fl.  =  AA,  AA,  AA,.  Every  cell 
can  accomodate  at  most  N,  electrons,  with  N,.  —  2U,  V/ 
Stt'.  where  the  factor  of  2  accounts  for  the  electron  spin. 
For  example,  if  we  take  n  =  lO''  cm  ’,  A  =  10'',  and 
AA,  =  AA,  =  AA  =  2  X  10'  cm  '  (A>  =  2.4.1  x  lO^cm  ' 
at  77  K).  we  would  have  F'  =  10  "  cm'  and  A,  =  6.45. 
A,  constitutes  the  maximum  occupancy  of  a  cell  in  the 


momentum  space  grid.  A  distribution  function  is  defined 
over  the  grid  in  momentum  space  by  counting  the  number 
of  electrons  in  each  cell.  The  distribution  function  so  de¬ 
fined  is  normalized  to  unity  by  dividing  the  number  in 
each  cell  by  A,,  for  use  in  the  rejection  technique.  We  note 
that  A.,  must  be  sufficiently  large  that  the  roundoff  to  an 
integer  does  not  create  a  significant  statistical  error.  In  the 
calculations  discussed  below,  we  used  AA.  =  8  x  Kf*  cm" ' 
and  AA^  =  AA,  =  8  x  lO',  which  assured  that  A,,  was 
greater  than  8.  Values  of  N  larger  than  10“*  were  tried  to 
check  the  convergence,  and  no  significant  differences  were 
found. 

When  a  .scattering  event  occurs  in  the  EMC  simulation, 
the  cell  <•,  containing  the  final  state  for  the  selected  tran¬ 
sition,  is  found.  The  normalized  distribution  function 
corresponding  to  that  cell,  is  compared  with  a  number  r 
randomly  chosen  between  0  and  1.  If  r  >  /.,  the  transition 
is  accepted  and  the  occupancy  of  cell  c  is  increased  by 
one.  If,  instead,  r  <  the  scattering  event  is  treated  as 
an  additional  self-scattering;  it  does  not  change  the  elec¬ 
tron  wave-vector.  This  corresponds  to  the  introduction  of 
an  auxiliary  self-scattering  mechanism  with  probability 
proportional  to  [  1  —  /(A')).  This  rejection  procedure  must 
be  included  since  this  latter  factor  was  ignored  in  the  ear¬ 
lier  steps  of  the  EMC  algorithm.  In  practice,  however,  one 
should  keep  the  cell  filling  factors  in  use  with  the  latter 
factor  in  calculating  the  scattering  rates,  but  the  results 
.should  be  statistically  the  same. 

The  relation  between  real  space  volume,  electron  con¬ 
centration,  and  the  number  of  simulated  electrons  guar¬ 
antees  the  correct  normalization  for  the  process.  This  se¬ 
lection  procedure  is  repeated  for  each  scattering  event. 
The  macroscopic  distribution  function,  that  is,  the  occu¬ 
pation  of  each  cell  of  the  grid,  is  updated  at  regular  inter¬ 
vals  (normally  of  the  order  of  the  typical  collision  time), 
and  the  procedure  is  iterated  until  the  stationary  distri¬ 
bution  is  obtained.  As  is  always  the  case  for  an  EMC  sim¬ 
ulation,  the  algorithm  just  described  allows  the  analysis  of 
the  transient  evolution  of  any  given  initial  distribution  to 
the  steady-state  one. 

In  the  next  section,  the  effect  of  the  Pauli  principle  on 
the  EMC  results  will  be  investigated.  Furthermore,  the 
influence  of  temperature,  electron  and  impurity  concen¬ 
trations,  and  electric  field  on  the  transport  properties  of 
degenerate  GaAs  will  be  analyzed.  The  ensemble  is  com¬ 
posed  of  10'*  electrons  and  only  the  steady-state  condition 
will  be  dealt  with.  All  of  the  results  will  be  for  EMC  sim¬ 
ulations  of  GaAs,  using  the  two-valley  model  described 
elsewhere  [10].  The  fields  used  are  sufficiently  low  that 
there  is  no  transfer  to  the  third  set  of  valleys.  The  L-val- 
leys  are  retained,  however,  because  there  is  a  slight  pos¬ 
sibility  that  the  electron-electron  interaction  can  kick  some 
electrons  to  a  sufficiently  high  energy  that  they  will  see 
these  valleys  |10].  The  figures  below,  however,  show  that 
these  electrons  are  very  few  in  number  and  that  there  is 
no  significant  electron  density  above  an  energy  of  0.1  cV, 
which  is  well  below  the  threshold  for  transfer  to  the 
L- valleys. 
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Fig.  1.  (a)  Energy  distribution  function,  (b)  energy  occupation  number,  and 
(c)  momentum  distribution  function,  for  n-GaAs,  with  (starred  curve) 
and  without  (solid  curve)  the  Pauli  exclusion  principle.  The  electric  field 
is  900  V/cm.  at  7”  =  77  K.  The  electron  (and  impurity)  concentration  is 
5  X  10”  cm  '. 

III.  Results 

Fig.  1  shows  the  effect  of  the  exclusion  principle  on  the 
electron  dynamics.  The  continuous  curves  refer  to  the 
classical  case;  that  is,  they  are  obtained  with  the  standard 
EMC  simulation.  The  starred  curves  indicate  the  results 
of  the  inclusion  of  degeneracy  in  the  manner  described  in 
the  previous  section.  An  electric  field  of  900  V/cm,  along 
the  z-axis,  is  used.  The  lattice  temperature  is  77  K,  and 
the  electron  and  impurity  concentrations  are  taken  to  be 
equal  to  5  x  lO'^  cm'\ 

The  energy  distribution  function /( £)  in  Fig.  1(a)  is  the 
number  of  electrons  in  an  energy  interval  AE  centered  at 
energy  E.  The  distribution  functions  here  are  normalized 
to  yield  the  integrated  value  n,  i.e.,  they  give  the  same 
weighted  density  at  E  and  are  plotted  in  arbitrary  units. 
This  integration  and  normalization  allows  one  very  clearly 
to  examine  the  shifts  of  particles  by  the  changes  in  height 


(or  area)  at  an  energy  value  (the  vertical  axis  is  arbitrary 
units— a. u.).  The  occupation  number  n(£)  is  shown  in  Fig. 
1(b);  this  is  obtained  by  normalizing /(£)  with  respect  to 
the  number  of  allowed  states  in  the  corresponding  energy 
cell,  and  relates  more  closely  to  the  quantity  /.  determined 
in  the  previous  section  although  we  are  plotting  the  results 
here  as  a  function  of  energy  rather  than  momentum.  The 
momentum  distribution  function  f(k,),  for  the  component 
of  the  electron  momentum  parallel  to  the  field,  is  the  pro¬ 
jection  of  f(lc)  onto  the  z-axis,  and  is  shown  in  Fig.  1(c). 
The  same  normalization  is  adopted  throughout  this  sec¬ 
tion  for  the  various  distribution  functions  in  order  to  allow 
comparison  between  different  figures. 

As  a  consequence  of  the  reduced  maximum  occupancy 
of  each  allowed  state,  electrons  are  pushed  into  the  tail  of 
the  distribution.  This  shift  corresponds  to  an  increase  of 
the  electron  “temperature”  of  about  50  percent  (from  140 
to  210  K)  when  the  Pauli  principle  is  included.  The  elec¬ 
tron  temperature  T,  (in  reality  the  average  random  kinetic 
energy  of  the  electron)  is  defined  in  the  usual  manner  for 
Monte  Carlo  procedures  as  that  part  of  the  total  kinetic 
energy  that  does  not  include  the  drift  component,  e.g., 

T,  =  {h^/Msm)  Ikl  +  kl  +  {k,  -  <k,)\  (2) 

This  corresponds  to  a  temperature  only  for  a  nondegener¬ 
ate  Maxwellian  distribution  function.  In  the  following,  all 
references  to  temperature  are  to  the  average  random  ki¬ 
netic  energy  of  (2),  and  should  be  so  interpreted.  (We  re¬ 
tain  this  usage  because  of  its  widespread  use  in  Monte 
Carlo  procedures,  but  admit  that  this  use  of  the  word 
“temperature”  is  not  an  accurate  descriptor.  In  fact,  this 
quantity  accounts  for  an  increase  of  the  electrochemical 
potential  due  to  degeneracy  plus  a  real-temperature  ef¬ 
fect.)  The  increase  of  the  electron  temperature  is  also  re¬ 
flected  in  the  broadening  of  the  momentum  distribution 
f(k.)  in  Fig.  1(c).  On  the  other  hand,  only  a  small  reduc¬ 
tion  of  the  drift  velocity  is  observed  in  the  degenerate  sit¬ 
uation.  This  is  due  to  the  increased  number  of  phonon 
emission  events,  as  more  electrons  are  pushed  above  the 
threshold  for  polar  optical  emission  processes  as  a  result 
of  the  exclusion  principle. 

At  lower  fields  (90  V/cm),  the  effect  of  degeneracy  is 
even  more  pronounced,  as  is  shown  in  Fig.  2.  The  elec¬ 
trons  have  a  very  small  drift  velocity  and  tend  to  remain 
in  the  lowest  portion  of  the  conduction  band,  with  an  av¬ 
erage  temperature  nearly  equal  to  the  lattice  temperature 
(continuous  curves)  in  the  nondegenerate  case.  Since  it  is 
impossible  to  allocate  more  than  two  electrons  per  state, 
a  migration  of  the  electrons  toward  higher  energy  states 
occurs  (starred  curves)  for  the  case  of  degeneracy.  In  the 
latter  case,  a  very  high  number  of  transitions  that  would 
normally  be  allowed  with  classical  statistics  are  instead 
forbidden  and  thus  rejected  due  to  the  occupation  of  the 
final  state.  In  the  simulations  illustrated  in  Fig.  2,  the  Pauli 
exclusion  principle  prevented  more  than  80  percent  of  the 
almost  lO'  scattering  events  simulated.  In  this  latter  case, 
the  electron  temperature  reaches  a  value  of  170  K,  as  com¬ 
pared  with  80  K  for  the  nondegenerate  EMC.  In  accor- 
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Fig  2  lul  Fncrgv  distrihuiion  function,  (b)  energy  iKTCupation  number,  and 
(c)  momenlum  disiribuiion  function,  for  n-GaAs.  with  (starred  curve) 
and  wiihoul  i solid  curve)  ihe  Pauli  exclusion  principle.  The  electric  field 
IS  90  V  cm.  dl  /■  -  77  K  The  eleclron  (and  impurity)  concentration  is 
.5  X  It)'  cm  ' 

dance  with  the  result.s  for  higher  helds,  the  drift  velocity 
shows  only  a  small  decrease.  In  Table  I,  we  give  the  values 
of  the  drift  velocity,  average  energy,  and  electron  '"tem¬ 
perature"  for  all  simulations  discussed  in  this  section. 

For  the  electric-held  strengths  considered  in  this  study, 
impurities  and  polar-optical  phonons  constitute  the  major 
scattering  priKesses.  A  simulation  that  keeps  the  .same 
electron  concentration,  but  greatly  reduces  the  number  of 
impurity  centers,  gives  drastically  different  results.  Fig. 
3  shows  the  results  for  an  impurity  concentration  of  only 
5  X  lO'*  cm  while  the  other  parameters  (including  the 
electron  concentration)  are  left  unchanged.  Physically  this 
can  be  achieved  by  exciting  electrons  into  the  conduction 
band  of  a  lightly  doped  sample  of  GaAs  by  laser  excita¬ 
tion.  by  modulation  doping,  or  by  injection  over  a  barrier 
into  a  lightly  doped  region,  as  in  a  planar-doped  barrier 
device  ill).  In  this  .situation,  electrons  are  weakly  scat¬ 
tered  by  the  impurities  and  their  drift  velocity  shows  a 


TABLE  1 


Non'degenerate  case  Results 


n. 

imp 

F 

^0 

<E> 

% 

(cm  ^) 

(V/cm) 

(K) 

(cm/ sec) 

(K) 

(K) 

5xlo’^ 

900 

77 

8.7xt0^ 

139 

128 

5x10^^ 

900 

77 

2.4x10^ 

190 

101 

90 

77 

1  .1x10^ 

81 

81 

Degenerate  case 

Results 

imp 

F 

"o 

n 

''d 

<E> 

T 

e 

(cm 

iV/cm) 

(K) 

(cm'^) 

(cm/ sec) 

(K) 

(K) 

5xlo'’ 

900 

77 

.x,o'^ 

b.lxlO*’ 

009 

204 

5xto'^ 

900 

77 

5X,0” 

2.2x10^ 

090 

.>1  S 

5x.o'' 

90 

77 

7.^X10 

175 

175 

5x,o” 

900 

JOO 

5X10” 

0 

4.1X10 

400 

396 

Ixio'^ 

900 

77 

Ixio'^ 

1  .5x10' 

155 

120 

very  sharp  increase.  This  shift  toward  higher  momenta  is 
evident  in  Fig.  3(c),  where  the  momentum  distribution 
functions  are  compared  with  those  for  the  heavily  doped 
case  (Fig.  1(c)).  As  a  consequence  of  the  higher  drift  ve¬ 
locity.  the  probability  of  finding  electrons  at  high  energies 
is  also  increased,  as  shown  in  Fig.  3(a)  and  (b).  Both  for 
degenerate  and  for  nondegenerate  statistics,  the  energy 
distribution  function  is  broadened.  The  effect  of  the  ex¬ 
clusion  principle  is  reduced  with  respect  to  the  heavily 
doped  case,  due  to  the  already  reduced  occupancy  of  the 
low-energy  states  by  the  hotter  and  faster  electrons. 

It  is  interesting  to  note  that  the  kink  present  in  the  con¬ 
tinuous  curve  (classical  case)  of  Fig.  3(a)  at  the  polar-op¬ 
tical  phonon  energy  (35  meV)  is  completely  washed  out  in 
the  degenerate  case.  This  kink  is  the  result  of  a  strong 
emission  of  optical  phonons  from  the  streaming  tail  of  the 
distribution,  which  creates  a  strong  push  of  electrons  by 
phonon  emission  back  into  the  low-energy  portion  of  the 
distribution.  The  streaming  character  of  this  process  is 
also  evident  in  the  elongated  tail  of  the  momentum  distri¬ 
bution  at  low  momenta  (continuous  curve  in  Fig.  3(c)). 
The  Pauli  principle  compensates  this  process  by  forcing 
the  electrons  back  into  the  high-energy  regions.  In  a  way, 
it  has  exactly  the  same  effect  on  the  energy  distribution 
function  as  the  electron-electron  interaction  IIOJ. 

Fig.  4  i.solates  the  effect  of  impurity  scattering  in  the 
degenerate  case.  The  continuous  and  the  starred  curves 
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Fig.  3.  (a)  Energy  di.\tribuii()n  function,  (b)  energy  occupation  number,  and 
(c»  momentum  distribution  function,  for  n-GaA.s,  with  (starred  curve) 
and  without  (solid  curve)  the  Pauli  exclusion  principle  The  electric  field 
is  900  V/cm.  at  T  =  77  K.  A  low-impurity  concentration  is  considered. 


Fig.  4.  (a)  Energy  distribution  function,  (b)  energy  occupation  number,  and 
(c)  momentum  distribution  function,  for  n-GaAs.  including  the  Pauli  ex¬ 
clusion  principle.  The  electric  field  is  900  V/cm,  at  T  =  77  K.  The  starred 
curve  and  the  solid  curve  are  for  high-  and  low-impurity  concentration, 
respectively. 


* 


I 


are  for  high  (5  x  lO'^  cm  ’)  and  low  (5  x  10'‘  cm  ’) 


concentration.s.  re.spectively.  The  drift  velocity  increases 
from  6  X  10*’ cm/s  to  2.2  x  10^  cm/s,  as  is  also  indicated 
by  the  shift  in  the  momentum  distribution.  The  increment 
in  the  drift  component  results  in  a  sharp  increase  in  the 
total  average  energy,  while  the  electron  temperature  is 
changed  only  slightly  (.see  also  Table  I).  This  velocity  in¬ 
crease  arises  primarily,  it  is  felt,  from  the  large  increase 
in  the  mobility  that  accompanies  the  decrease  in  doping 
level.  This  difference  is  significant,  as  it  is  now  clear  that 
attempts  to  study  the  electron-electron  interaction  by  using 
heavily  doped  .samples  can  fail  due  to  the  masking  effects 
of  the  impurity  scattering. 

The  comparison  of  the  (Kcupation  number  obtained  from 
the  present  calculation  and  the  equilibrium  Fermi  distri¬ 
bution  is  shown  in  Fig.  5.  The  continuous  curve  is  for  the 
equilibrium  theoretical  distribution,  while  the  starred  one 
is  for  the  degenerate  EMC  result.  The  Fermi  level  in  a 


>51  CJi'  i 


pation  number,  lor  n  GaAs,  at  77  K  The  elcetron  eoneeniralion  is  5 
10'  em  '.  and  the  eleelrie  field  9(K)  V/em. 


semiconductor  at  the  lattic  temperature  T„  is  found  from 

(3) 

where 


n  —  n,  F 1 2  lEf  (7',,)/kfl7',,] 
n.  =  0.25{2mksTJith')^' 


(4) 


. 
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Fig.  6.  Comparison  of  the  Fermi-Dirac  distribution  with  the  EMC  occu¬ 
pation  number,  for  n-GaAs,  at  300  K.  The  electron  concentration  is  S  x 
lO”  cm  ’,  and  the  electric  field  900  V/cm. 


Fig.  7.  (a)  Energy  distribution  function,  (b)  energy  (K'cupalion  number,  and 
(c)  momentum  distribution  function,  for  n-GaAs.  including  the  Pauli  ex¬ 
clusion  principle.  The  electric  field  is  900  V/cm.  the  electron  concentra¬ 
tion  .5  X  lO”  cm  '.  The  starred  curve  is  for  T  =  3(K)  K;  the  solid  one 
IS  for  /•  =  77  K. 


is  the  effective  density  of  states  in  the  conduction  band, 
and  F,  2  is  the  Fermi-Dirac  integral.  In  fact,  the  random 
energy  given  by  T,  should  be  used  in  (3)  and  (4).  and  this 
further  reduction  in  the  Fermi  energy  is  apparent  in  Fig. 
5.  The  EMC  curve  is  the  same  as  shown  in  Fig.  1(b).  At 
77  K,  the  chemical  potential  is  27  meV  (at  T,  =  200  K) 
above  the  bottom  of  the  conduction  band.  The  noise  in  the 
lowest  part  of  the  EMC  curve  derives  from  the  small  num¬ 
ber  of  states  enclosed  in  the  first  energy  cells,  due  to  the 
discretization  of  ?-space  described  before.  The  higher 
population  at  high  energies  of  the  EMC  curve,  relative  to 
the  equilibrium  curve,  corresponds  to  the  heating  of  the 
electrons,  whose  temperature  is  almost  three  times  higher 
than  the  lattice  temperature  under  these  conditions. 


At  increasing  ambient  (or  electron)  temperatures,  the 
chemical  potential  moves  toward  the  bottom  of  the  con¬ 
duction  band.  The  occupation  probability  of  the  low-en¬ 
ergy  states  decreases.  Fig.  6  gives  the  equilibrium  dis¬ 
tribution  (continuous  curve)  and  the  EMC -derived 
distribution  (starred  curve)  at  room  temperature.  In  this 
case,  Ef  is  equal  to  I.I  meV.  The  effect  of  the  higher  ther¬ 
mal  motion  of  the  electrons  at  room  temperature  is  also 
evident  by  looking  at  the  results  shown  in  Fig.  7,  where 
the  degenerate  EMC  distribution  function  at  300  K  (con¬ 
tinuous  curve)  is  compared  with  that  at  77  K  (starred 
curve). 

As  the  electron  concentration  is  reduced,  the  classical 
situation  is  approached  so  that  the  exclusion  principle  no 
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